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Abstract. 

We demonstrate a stable, mobile, dipolar or nondipolar three-dimensional matter- 
wave soliton in the vortex core of a uniform nondipolar condensate. All intra- and 
inter-species contact interactions can be repulsive for a strongly dipolar soliton. For 
a weakly dipolar or nondipolar soliton, the intra-species contact interaction in the 
soliton should be attractive for the formation of a compact soliton. The soliton can 
propagate with a constant velocity along the vortex core without any deformation. 
Two such solitons undergo a quasi-elastic collision at medium velocities. We illustrate 
the findings using realistic interactions in a mean-field model of binary 87 Rb- 85 Rb and 
87 Rb- 164 Dy systems. 
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A bright soliton is a self-bound object that travels at a constant velocity without 
deformation in one dimension (ID), due to a cancellation of nonlinear attraction and 
defocusing forces. An ID dark soliton is a dip in uniform density, which also moves with 
a constant velocity maintaining its shape. ID Solitons have been observed in nonlinear 
optics pp, and in Bose-Einstein condensate (BEC) [2]. Experimentally, bright matter- 
wave solitons were created in BECs of 'Li [3j and 85 Rb atoms [Tj. Dark solitons were 
also observed in BECs of 8 'Rb [5] and 23 Na [6]. The ID set up is obtained by putting 
confining traps in directions perpendicular to the motion of the soliton. However, 
a three-dimensional (3D) trap-less soliton cannot be formed for a cubic nonlinearity, 
generally encountered in BEC and nonlinear optics, due to collapse [7]. The collapse 
can be stopped in a weaker saturable [8], cubic-quintic PJ or quadratic [10] nonlinearity, 
or by an application of nonlinearity and/or dispersion management [11] . In nonlinear 
optics, ID temporal solitons [T2j as well as lattice solitons in arrays of nonlinear optical 
wave guide in ID [ 13] and in two dimensions (2D) [14] and 3D na, with mo diked 
dynamics/nonlinearity, have been observed. 

Here we demonstrate the formation of a trap-less matter-wave soliton in the core 
of a quantized vortex of a uniform nondipolar BEC, which we call a binary nondipolar 
vortex-soliton, which is a 3D analogue of an ID dark-bright soliton pL] studied previously. 
These solitons are shown to be stable and execute steady oscillation for very long time 
under a small perturbation. In the case of a 3D vortex-soliton all interactions can be 
repulsive except the intra-species interaction in the soliton. The soliton can swim freely 
with a constant velocity along the vortex core. Because of the strong localization of 
the soliton due to inter-species contact repulsion, the soliton can move without visible 
deformation, The collision between two integrable ID solitons is truly elastic DEI- 
However, at medium velocities the collision between two solitons is found to be quasi¬ 
elastic without visible deformation. As no modification of the nonlinear interactions is 
suggested such a trap-less soliton can be realized in a laboratory. 

The observation of dipolar BECs of 164 Dy [16], 168 Er [l7j and 52 Cr [18] has initiated 
studies of new types of BEC solitons. For example, one can have a dipolar BEC soliton 
for a repulsive contact interaction [T9] , in 2D [20] or on optical-lattice potentials [23]. 
These new dipolar solitons were possible due to the peculiar nature of dipolar interaction. 
The dipolar BEC solitons of a large number of atoms stabilized by a long-range dipolar 
attraction, could be robust and less vulnerable to collapse in the presence of a short- 
range contact repulsion [I9U22]. Hence, we also consider a dipolar soliton in a nondipolar 
vortex core, which we call a dipolar vortex-soliton. Actually, a strong short-range inter¬ 
species contact repulsion between the atoms of the matter-wave soliton and the atoms 
of the vortex localizes the soliton in 3D. 

In the present investigation we use a mean-held model described in section [2] First 
we present an ID model in section 2.1 useful for an analytic understanding of the 


formation of a dark-bright soliton. The 3D mean-held model is presented in section 2.2 
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In section [3] we present numerical results for the formation of a binary vortex-soliton. 
A summary of our findings is given in section [4j 


2. Mean-field model 


2.1. One Dimension 


A vortex in a uniform BEC bears similarity with an ID dark soliton in having a 
hole along the axial z direction and is often called a 3D dark soliton [I]. A binary 
vortex-soliton is the 3D analogue of the well-known ID dark-bright soliton. Hence, to 
understand how a 3D vortex-soliton can appear, we consider the following integrable 
binary ID dark-bright soliton model of atoms of same mass m, same inter- and intra¬ 
species scattering lengths a and same number of atoms N 

aj - (x,t) r Ti 2 d 2 2 


ih- 


dt 


2 m dx 2 


+ gJ2 l&l 2 


2—1 


( 1 ) 


where i,j = 1,2 , g — 2ah 2 N/(md 2 ), d p = \JhJ(nud), where u is the frequency of a 
strong harmonic trap in the binary mixture in the transverse directions. Scaling the 
wave functions by |0j| 2 = g\(fi\ 1 we obtain 


.dfjj(x, t ) 
' dt 


1 8 2 

2 dx 2 




( 2 ) 


in units h = m = 1. Equation ([2]) with all-repulsive interactions has the analytic dark- 
bright soliton [lj 


^i(x, t) = /3tanh [a(x - vt)]e ivx ~^ 2/2+/?2)t , (3) 

0 2 (z, t) = ysech [a{x - vt )\ e ivx+i tt a2 -v 2 )D-P 2 }f (4) 


where a and (3 (/3 > a) are constants which control the intensity and width of the 
solitons, 7 = \Jfi 2 — a 2 and v is the velocity. The appearance of bright soliton Q in 
the all-repulsive equation (J2| is counterintuitive and is possible due to the coupled dark 
soliton ([3]). Without losing generality we impose the normalization / \if 2 (x,t)\ 2 dx = 1, 
yielding (3 = \Ja 2 + a/2. For v = 0 the solutions have the form 

V’i(^,t) = \joL 2 + a/2tanh(crr) exp[— i(a 2 + a/2)t], (5) 

-02 {x,t) = \/ a / 2sech(ax) exp [—i(a 2 + a)t/2]. (6) 

To solve ([2]) numerically the simulation was performed in a ID box in the domain 
x = ±50. In figure [l] we plot the numerical matter-wave densities of the dark-bright 
soliton and the analytic results ([5]) and (j6|.The bright soliton sits in the central hole 
of the dark soliton and the inter-species repulsion between the (outer) dark and (inner) 
bright solitons confines the latter. Similarly, the soliton of a vortex-soliton can be 
confined in the radial x — y plane by the inter-species repulsion between the vortex 
and soliton. The conhnement along the z direction is obtained by the intra-species 
dipolar attraction with the dipoles polarized along z axis and/or by intra-species contact 
attraction. 
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Figure 1. (Color online) Densities Pj(x) = \i/jj(x)\ 2 from analytic (ana) model |5j) and 
([HJ) with a = 0.198, /? = y/a 2 + a/2 = .372 compared with those of numerical (num) 
simulation by imaginary-time propagation in the space domain x = ±L = ±50. The 
normalization of the wave functions are J __ L pi(x)dx = 12.42 and f x __ L p 2 {x)dx = 1. 


2.2. Three Dimensions 


We will consider a matter-wave vortex-soliton in the form of a vector soliton without 
any trapping potential and present the binary BEC model appropriate for this study. 
The first component (j = 1), with the vortex, is nondipolar and the second component 
(j = 2), with the soliton, can be dipolar or nondipolar. As in the ID dark-bright soliton, 
the intra-species repulsion prevents the matter-wave soliton from escaping radially. The 
dipolar interaction or intra-species contact attraction prevents the soliton from escaping 
in the axial direction. The mass, number of atoms, and scattering length for the two 
species are rrij,Nj,aj, respectively. The intra- (Vf) and inter-species (V 12 ) interactions 
for atoms at r and r' are [24J 


= 4^a^(R) (R) = 2rfq 12 yR) 

m 1 ^ m R 

3a dd rr 2 D dd (R) 4Trti 2 a 2 5(R) 

— -1-, 

m-2 m-2 


(7) 

( 8 ) 


where R = (r — r'), a dd = poAi 2 ' m 2/(127r^. 2 ), Vd d (R) = (1 — 3cos 2 6 I )/R 3 , the reduced 
mass m R = m\m 2 /{m\ + mf), a± 2 is the inter-species scattering length, a dd is a dipolar 
length to measure the strength of dipolar interaction, /i 0 is the permeability of free 
space, n is the magnetic moment of each atom in the dipolar soliton, 9 is the angle 
made by the vector R with the polarization z direction. The dimensionless mean-field 
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Gross-Pitaevskii (GP) equations for the trap-less binary mixture are [21] 


■ <90i(r, t) 

1 dt 

■ d<j> 2 (r,t) 

1 dt 

+ gdd [ V r d d(R)|02(r , ,t)| J (ir , |02(r,t), 


- Y + #l|0l| 2 + ^712 1021 2 

V 2 

— m l2^Z -h ^21021 2 + <7211011 2 


(9) 


( 10 ) 


where mi 2 = mi/m 2 , r/i = 47raiiVi, g 2 = 47ra 2 lV 2 'mi2, ^i 2 = 27rmia 12 A r 2 /m fi , 
r/ 21 = 27rmia 12 iV 1 /m R , g dd = 31V 2 r? 'T2 ( hM/m^. In (J9| and (10), length is expressed in 
units of a scale l , probability density |0j| 2 in units of C 3 , energy in units of h 2 / {2m\l 2 ) 
and time in units of t 0 = 2 mil 2 /h. 

To hnd a stationary quantized vortex of angular momentum £ along z axis in 
component 1, also called a dark soliton with circular symmetry, we look for axially- 
symmetric solution 4> 1 (r,t) in x — y plane: 0i(r ,t) = $i(r ,t)e lClf , where ip is the 
azimuthal angle and $i(r,f) satisfies [I] 

.9$i(r ,t) 


dt 


1 / d 2 d 2 d 2 

2 ' dz 2 dx 2 + dy 2 


+ C 2 /(x 2 + y 2 ) + g 1 \$ 1 \ 2 + gn\fa\ 2 }$ 1 (r 1 t), (11) 

with the boundary conditions <l>i(x = 0 ,y = 0 ,z) = 0,&i(x —> oo,y — > oo,^;) = 
constant. The same on the ID dark soliton (J3| are very similar: ipi(x = 0) = — » 

oo) = constant. For a matter-wave soliton of component 2 in the vortex core of 
component 1 we solve the axially-symmetric equations © and ( floj ) . We consider 
a vortex of unit circulation £ = 1. 


3. Numerical Results in three dimensions 

In the 3D simulation of the binary vortex-soliton we consider the nondipolar 8 'Rb- 85 Rb 
and the dipolar 87 Rb- 164 Dy mixtures. The 164 Dy atom has the magnetic moment /i = 
10 ys [16] with p B the Bohr magneton so that the dipolar length cidd( 164 Dy) ~ 132.7a 0 
with a 0 the Bohr radius. We use scattering lengths a( 8 'Rb)= a( 8 'Rb- 85 Rb) = a( 8 'Rb- 
164 Dy) = 120ao and take a( 85 Rb) and a( 161 Dy) as variables. The experimental values of 
these scattering lengths are not known precisely. The exact values of the inter-species 
scattering lengths are not important for our analysis. These positive scattering lengths 
are used as they simulate the inter-species repulsion required for the formation of binary 
vortex-soliton. Furthermore, if needed, the variation of the scattering lengths can be 
achieved by the Feshbach resonance technique [25]. We solve (10) and © by the split- 
time-step Crank-Nicolson method using both real- and imaginary-time propagations in 
Cartesian coordinates using a space step of 0.2 ~ 0.4 and a time step of 0.0025 ~ 0.005 
[2>J. The dipolar term is treated by a Fourier transformation in momentum space using 
a convolution theorem In all cases we take the length scale l = 1 /irn and time 
scale t 0 = 2m( 8 'Rb )/ 2 /% = 2.74 ms. The numerical simulation is performed in a cubic 
box, limited by x = y = z = ±50, containing 400000 8 'Rb atoms of density 4 x 10 11 /cc. 
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Figure 2. (Color online) (a) Integrated ID density P\d{x) and piu(z) of a matter- 
wave soliton of 1000 S5 Rb atoms in the vortex core of a uniform 87 Rb BEC of 
4 x 10 5 atoms in a box of volume 100 3 /im 3 . The same of (b) 500, (c) 1000, and 
(d) 20 00 164 Dy atoms, (e) 3D density P3d{x,0, 0) of the vortex core for (a) and 
(d) together with analytic fit (121 and healing length estimate (healing) p(x, 0,0) = 
0.00000103a; 2 /(a: 2 + 2£ 2 ) with healing length £ ss 4.072 pm. (f) Integrated 2D density 
P 2 d(x, y ) of the 87 Rb vortex and the 164 Dy soliton of (d). The parameters are a( 8 'Rb- 
85 Rb) = a(Rb-Dy) =120ao- Variables in all figures are dimensionless. 


We consider a matter-wave soliton of 1000 nondipolar 85 Rb atoms in the 8 'Rb 
vortex core and perform imaginary-time simulation. A large intra-species attraction 
with scattering length a( 85 Rb) = —80a 0 was necessary to obtain a compact soliton of 
small size. In the nondipolar case, unlike in ID equation (J2|, no 3D vortex-soliton can be 
obtained for repulsive inter- and intra-species interactions. An attractive intra-species 
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interaction facilitates the formation of the soliton. For smaller values of intra-species 
attraction, the size of the soliton was larger resulting in computational difficulty. The 
ID densities of the nondipolar soliton, defined by Pid{x) — 5 dy f dz |£>(r)| 2 , etc., are 
plotted in figure [2^a). 

In the dipolar case, we consider solitons of 500, 1000, and 2000 164 Dy atoms, in 
the 8 'Rb vortex core, with intra-species scattering lengths a( 164 Dy) = —40a 0 ,10a 0 , and 
30ao, respectively. Again the exact values of the scattering lengths are not important. 
These three values of scattering lengths simulate three distinct interactions: attractive, 
weakly repulsive and moderately repulsive. If needed, the intra-species and inter¬ 
species interactions can be manipulated by independent optical and magmetic Feshbach 
resonances in a laboratory ra- The corresponding ID densities of these solitons are 
presented in Figs. [ 2 ] (b), (c), and (d), respectively. The 3D density p 3 £>(r) = | < h 1 (r,t)|- 
of the axially-symmetric vortex core is plotted in figure [2](e) . We also show the 
variational healing-length estimate of this density for an isolated vortex without the 
soliton: p(x,y, 0) = po(x 2 + y 2 )/(f 2 + x 2 + y 2 ), where po is the density away from the 
vortex core, and the healing length £ = l/^8npoa ~ 4.072 pm. The numerical solution 
of the GP equation for the isolated vortex agrees well with healing-length estimate. The 
deviation of the vortex core of the binary vortex-soliton from the healing-length estimate 
of an isolated vortex is due to the presence of the soliton which increases the size of 
the vortex core by inter-species repulsion. In figure [2|f) the integrated 2D densities 
p 2 D(x,y) = f dz\(f(r)\ 2 for the vortex and the soliton corresponding to figure |2^d) are 
plotted. Figures |2](a) — (f) show that the soliton is localized in the vortex core given by 
the minimum of the vortex density plotted in figure |2](e) . 

A qualitative understanding of the formation of the matter-wave soliton can be 
obtained if we note that the densities corresponding to the vortex wave functions as 
presented in figure [2|e) for the four cases studied above are practically the same given 
by the following function 

<P?(r) = A(1 - e~ ix2+y2)s )(l + u- ue~ z2s ), (12) 


with A = 0.00000103,(5 = 0.01, v = 0.25. The function (12) with v = 0 simulates a 
vortex in a uniform BEC. A non-zero value of v includes the effect of a small distortion of 


the vortex inside the soliton due to the presence of the soliton. The analytic density (12) 
is also plotted in figure [2|e) in good agreement with the numerical result. We perform 
an approximate variational calculation for the formation of a soliton with vortex density 
(12) substituted in (10) using the following Gaussian trial function for the soliton: 

<f>2 = 


7 r 


-3/4 


Wnx/W 


exp 


x 2 + y 2 




z 

2wV ’ 


(13) 


where w p and w z are the widths. The Lagrangian (10) can be written as 




dr 


mi2A r 2 | V£>2|" + N2gif>\ + 2A r 2 p2i| < hi| 2 |02 


dr / dr , A r 2PddDdd(R.)|02(r)| J |0 2 (r , )P 


( 14 ) 
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Figure 3. (Color online) 3D isodensity contour of the binary vortex-soliton, showing 
the vortex core (grey, pink in color) and the soliton (black, blue in color) profiles, 
corresponding to (a) figure [2](a) , (b) figure J2|b), (c) figure [2|c), and (d) figure [2](d). 
Density on the contour is 10 b /cc. 


N 2 mu r 1 1 


L W 


+ 


N 2 g 2l ASw 


e 2 (w. 


2 w 2 


1 + u 


+ 


Nfm 12 [a 2 - a dd f(K,)] 


\[2fnw 2 w z 


e(w z ) i 


(15) 


where k = w p /w z , /(«) = [1 + 2 k 2 —3 n 2 dfn)]/ (1 —k 2 ), d(n) = arctan(\/ k 2 — V)/\J k 2 — 1 , 
e(x) = \Jl + 8x 2 . Minimizing Lagrangian (15) we find the following conditions to 
determine the widths of the soliton mm 


1 : N 2 e(n ) 2w p g 2 iA5 - 

e A (w p )m l2 L 


+ 


p 

77*12 

wl 


1+1/- 


+ 


\[2fnw 3 w z 


2N 2 m l2 h{n) 2g 21 Aisd 2 w 2 w z 
\/2nw 2 w z 


2 w 2 
7k 2 


e 2 (w p )e 3 (w z 


e(w z ) 
= 0 , 


= 0 , 


(16) 


(17) 


where e(«) = 2a 2 — add[2 — 7k 2 — 4k 4 + 9k 4 <7(k)]/(1 — k 2 ) 2 , h(n) = a 2 — add[l + 10 k 2 — 
2k 4 — 9k 2 <7(k)]/(1 — k 2 ) 2 . A solution of ( jl 6 | ) and (17) determines the widths, and hence 
the sizes of the soliton. For the solitons of figure [2] (a) — (d) the numerical and analytic 
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root-mean-square (RMS) sizes are given in table [lj The analytic result presented is not 
a variational solution of the full dynamics given by GD and ( [To] ) and hence there are no 
variational bounds on the energy or sizes. Nevertheless, the approximate analytic results 
of table [lj in reasonable agreement with the numerical results, provide a qualitative 
understanding of the formation of the soliton. 

The density for the dipolar matter-wave solitons is strongly anisotropic with distinct 
ID densities Pid(x) and Pid(z) as can be found in Figs. |2](b) — (d), whereas in the 
nondipolar case, shown in figure |2](a), these two densities are nearly equal. The 
anisotropy in the shape of the soliton arises partly due to the anisotropy of the dipolar 
interaction and partly due to the anisotropy of the vortex core. The increase in total 
dipolar interaction makes the soliton more elongated (prolate) along the polarization z 
direction. This is reflected in the ratio of numerical RMS sizes (z)/(x) shown in table 

m 


The 3D isodensity contours of the four matter-wave solitons of Figs. ga)-(d) 
together with the respective vortex cores are shown in Figs. [3] (a) —(d), respectively. 
The central prolate spheroid is the soliton and the outer shell is the vortex core. The 
net dipolar interaction and hence the anisotropy of the soliton increases with the number 
of dipolar 164 Dy atoms. The 3D profile of the soliton is more prolate, in Figs. |3](b) — (d), 
due to dipolar interaction, compared to that in figure |3](a) for the nondipolar case. 
The nondipolar soliton is also slightly prolate due to the axially-symmetric inter-species 
repulsion of the vortex. 

By carefully adjusting the parameters of the mean-held model equations a balance 
between attraction and repulsion can be achieved to obtain a vortex-soliton similar 
to a Townes soliton in 2D [2Sj, which is weakly unstable and collapses upon a small 
perturbation [7]. So it is quite relevant to establish the stability of the vortex-soliton. 
To demonstrate the stability of the soliton, we consider the one in figure [3|a) and subject 
the stationary state(s) obtained by imaginary-time propagation to real-time propagation 
introducing a small perturbation, e.g., jumping intra-species scattering length a( 85 Rb) 
from —80a 0 to —80.5ao at t — 0. Long-time stable oscillation of the resultant RMS 
sizes, illustrated in figure [4] (a), guarantees the stability of the soliton. To test the 
stability of the soliton after a small displacement along x direction, we performed real¬ 
time simulation of the vortex-soliton of figure [3] (a) after displacing it through 2 units of 
length along x direction. The soliton is found to come back to its stable position on the 


Table 1 . Numerical and analytic RMS sizes of the four solitons presented in Figs. 


|3N— (d). 



ll')num 

5 y)anal 

ip) num 

\Z) anal 

(z / X s ) num 

(a) 

4.89 

4.762 

6.97 

7.519 

1.43 

(b) 

4.66 

4.341 

11.05 

10.557 

2.37 

(c) 

4.47 

4.226 

10.90 

11.711 

2.44 

(d) 

3.79 

3.870 

9.72 

11.884 

2.56 
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Figure 4. (Color online) (a) RMS sizes (x),(z) during breathing oscillation of the 
nondipolar soliton of Figs. [2ja) and^a) initiated by a sudden change of a( 85 Rb) from 
—80a o to —80.5ao- (b) Initial (i) and final (f) ID densities along x and z directions 
after real-time dynamics of the vortex-soliton of figure [3] (a) during 1000 units of time. 
The dynamics is started after shifting the initial soliton along x direction through 2 
units of length and maintaining the initial vortex position unchanged. In both (a) and 
(b) the initial state was the imaginary-time stationary profile of the vortex-soliton of 
figure^ a). 


z axis at the center of the vortex core without executing oscillatory motion along x axis. 
In figure [4] (b) the plot of initial and final ID densities of the soliton after 1000 units of 
time confirms the stability. This transverse stability of the vortex-soliton is important 
for an experimental realization. 

Next we study the head-on collision between two solitons moving along the vortex 
core in opposite directions. The imaginary-time profile of the vortex-soliton of figure 
|3][a) is used as the initial function in the real-time simulation of collision, with two 
identical solitons placed at z — ±50. To set the solitons in motion along the z axis 
in opposite directions the soliton wave functions are multiplied by exp (±ivz),v = 20. 
To illustrate the dynamics upon real-time simulation, we plot the time evolution of ID 
density pio(z,t) in figure [5j(a) and its contour plot in figure |5](b). The same for the 
collision of two solitons of figure [3j(c) are shown in Figs. [5](c) and (d). In Figs. [5j(a) 
and (b), the dimensionless velocity of a soliton is ~ 2.5, whereas in Figs. |5](c) and (d) 
this velocity is ~ 1.25. The quasi-elastic nature of collision is established at a relative 
velocity (two times the velocity of a single soliton) of about 2 ~ 5 in dimensionless units. 

4. Summary 

Summarizing, we demonstrated the possibility of the creation of a nondipolar or dipolar 
matter-wave soliton in the vortex core of a uniform nondipolar BEC. The soliton is 
localized by a strong inter-species repulsion. This binary vortex-soliton is a stable 
stationary state. A dipolar soliton can be created for repulsive inter- and intra-species 
contact interactions. However, for the creation of a nondipolar soliton an attractive 
intra-species contact interaction in the soliton is necessary. The matter-wave soliton 
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Figure 5. (Color online) (a) The ID density pio(z,t) and (b) its contour plot during 
collision of two nondipolar solitons of 100 0 85 Rb atoms of figure [ 3 ] (a) initially placed 
at 2 = ±50, upon real-time propagation. The initial wave functions are multiplied by 
exp(±*20z) to set them in motion. The same for two dipolar solitons of 1000 lb4 Dy 
atoms of figure [3] (c) are shown in (c) and (d). 


can move with a constant velocity along the vortex core without any deformation. The 
stability of the vortex-soliton is demonstrated by a stable oscillation of the soliton upon a 
small perturbation in real-time simulation using the initial states obtained in imaginary- 
time propagation, viz. figure [4j At medium velocities, the collision between the two 
solitons is quasi elastic with no visible deformation, viz. figure [5] 

The techniques of generating a vortex in BECs [29] are well known, hence a binary 
vortex-soliton can be realized in experiments. Here, for the sake of computational 
simplicity, we considered the vortex in a uniform BEC. However, for an experimental 
observation, a small matter-wave soliton can be created in the core of a large low-density 
weakly-trapped BEC vortex. To achieve this the binary vortex-soliton should be realized 
under a weak harmonic trap on both components in a laboratory and eventually the weak 
trap on the soliton should be slowly removed to obtain the trapless soliton in the vortex 
core maintaining the weak trap on the vortex. The weak initial harmonic trap on both 
components will bring the soliton inside the vortex core. We tested the stability of the 
binary vortex-soliton under a small transverse perturbation. The transverse confining 
force on the soliton due to inter-species repulsion brings the soliton back to the center 
inside the vortex core. An interesting future work would be to study the possibility of 
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creating a binary vortex-soliton in two components of a spin-orbit coupled BEC with 
one component supporting a coreless vortex [30] and the other supporting a soliton. 
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